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Abstract 

The goal of the SIMBA cohaboration is to provide a global fit to the available data in inclusive 
B — )• Xg^ and B — )• Xuiu decays. By performing a global fit one is able to simultaneously 
determine the relevant normalizations, i.e. the total B — )• rate and the CKM-matrix element 
\Vub\, together with the required input parameters, most importantly the 6-quark mass and the 
6-quark distribution function in the i?-meson, called the shape function. This strategy is analogous 
to the determination of \Vcb\ from global fits to inclusive B — >■ Xcii' decays. In this talk, we present 
preliminary results for the shape function and IC^'^'Vtb^sl' which parametrizes the total B — >■ Xgj 
rate, from a global fit to the available B — > Xsj measurements from BaBar and Belle. 
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1. A GLOBAL FIT APPROACH TO B ^ X^-y AND B X^iv 



The flavor- changing neutral current B — Xg'^ process is very sensitive to contributions 
from new physics beyond the Standard Model (SM). To exploit this sensitivity and constrain 
new physics, an accurate extraction of the B — )■ Xg'-y rate from data is desirable. Currently, 
the experimentally measured partial branching fractions with a cut on the photon energy, 

> E^^^, are extrapolated to the partial branching fraction for a fixed cut E^ > 1.6 GeV, 
which yields Q B{E^ > 1.6) = (3.55±0.24±0.09) x 10"^ This value is then compared to the 
fixed next-to-next-to-leading order (NNLO) SM prediction from Refs. j2|, 0], B{Ery > 1.6) = 
(3.15 ± 0.23) X 10~^. This procedure is used to compare experiment and theory because 
fixed-order perturbation theory can only be applied for low enough values of i?™*. However, 
since the extrapolation down to E'™* = 1.6 GeV uses theory, it still requires a theoretical 
calculation of the decay rate with the actual value of E^'^'^ used in each measurement. In 
other words, irrespectively of how the extrapolation is performed, the comparison of theory 
and experiment always happens effectively at the measured E™^. 

The currently performed extrapolation and average fi\ assumes a model for the leading 
shape function, which introduces a systematic uncertainty from the model dependence. The 
extrapolation also has an additional theory uncertainty which is correlated to some extent 
with that of the fixed-order SM prediction at i?^^* = 1.6 GeV. To minimize the effect of the 
extrapolation, from each experimental analysis typically the measurement with the smallest 
possible i?™* is used, which has the largest experimental systematic uncertainty. Hence, 
a small subset of the experimental information dominates the result, and in particular the 
more precise measurements at higher values of E^ cannot be utilized. 

Our strategy, which was proposed in Ref. avoids these drawbacks. Performing a 
global fit allows one to minimize the uncertainties by making maximal use of all available 
data at any E^. At the same time it allows for a consistent treatment of correlated un- 
certainties, both experimental and theoretical, as well as from input parameters. We use a 
model-independent treatment of the shape function, such that its shape and uncertainty is 
determined by the shape and uncertainties in the measured E^ spectra. The overall b ^ s'j 
transition rate, which holds the sensitivity to new physics, is parametrized by the combi- 
nation ICy^'^^V^bV^*!, defined below, and is determined by the normalization of the measured 
spectra. Its value obtained from the global fit can then be compared to its SM prediction 
in order to constrain possible contributions from new physics. In this way, measurements 
at all E^ contribute optimally to constrain the total B — )■ Xgj rate. Furthermore, the 
extracted shape function provides a necessary input for the determination of \Vub\ from 
inclusive B — )■ Xuiu decays, e.g. via a combined global fit, which is left for future work. 

2. THEORY 

A. Treatment of the Shape Function 

The shape function renormalized in MS, S{uj,n), which enters the description oi B ^ 



Xs7, can be factorized as 0] 




(1) 
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Here Co{u, fi) is the MS-renormahzed 6-quark matrix element of the shape- function operator 
calculated in perturbation theory, while F{k) is the nonperturbative contribution to S{u, /i). 
The hats on F{k), Co{uj, /i), and mj, below indicate that they are defined in a short-distance 
scheme. Here we use the IS" scheme, rhb = fnl^ , see Ref. [3] for more details. 

The construction in Eq. ([T]) has several advantages. It ensures that S{u, /x) has the correct 
perturbative tail at large u and also the correct /x dependence and RGE, which both come 
from Cq{(jJ,^). For small w, the shape of S{uj,fi) is determined by F{k). Hence, F{k) is the 
nonperturbative parameter that determines the shape of the B — Xgj spectrum at large 
E^ and which we need to extract from the data. For simplicity, we will refer to F{k) as the 
shape function in the following. In contrast to S{ui,jj,), F{k) falls off exponentially at large 
k, so the moments of F{k) exist without a cutoff, and information about rrib and matrix 
elements of local operators can be incorporated via constraints on the moments of F{k). For 
example, 

jdkF{k) = l, jdkkF{k) = mB-mb. (2) 
To fit F{k) from data we follow Ref. [4i] and expand it in a complete orthonormal basis, 

- CO , -1 2 /. oo 

^c„/„(-) With UkF{k) = Y,cl = l. (3) 



A 



1 - ' 2 

F{k) = - 



l-n=0 



n=0 



The basis functions fn{x) are given in Ref. ^] and A ^ Aqcd is a dimension-one parameter 
of the basis. Since the functional basis is complete, Eq. (jS]) provides a model-independent 
description of F{k), where its shape is parametrized by the basis coefficients c„. 

By fitting the coefficients c„ from data, the experimental uncertainties and correlations 
in the measured spectra are captured in the uncertainties and correlations of the fitted c^. 
In practice, the data only allow a fit to a finite number of coefficients, so the expansion 
must be truncated after -|- 1 terms. This introduces a residual model dependence from the 
chosen functional basis, in particular the value used for A. The overall size of this truncation 
uncertainty scales as 1 — "^^^q c\. The optimal values for A and are determined from the 
data. The value for A is chosen such that the fitted series converges quickly, and the number 

of fit coefficients should be large enough such that the truncation uncertainty is small 
compared to the experimental uncertainties of the fit coefficients. In other words, we let the 
available data determine the precision to which the functional form of the shape function is 
known, by including as many basis coefficients in the fit as possible given the available data. 
Hence, our approach allows for an experimental determination of the shape function which 
is model independent and yields reliable, data-driven uncertainties. 



B. Master Formula for B Xs'y 

The B — )■ photon energy spectrum is given by 



dr Gla _ , 

' ""b\VtbVts 



dE^ 27r 



F"em p3 --2 it/ t/*|2 



X I \C'rf I dA; W-jj{k)F{mB -2E^-k) + Y^ Wj-j,^ FmirriB - 2E^ 
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+ fdkY^ [2Re{C'^''^)Q Wnik) + QCj Wij{k) F{mB - 2E^ - k)\ . (4) 

The expressions entering Eq. will be given in Ref. jsj. The function Wn{k) contains 
the perturbative corrections to the b S'j decay via the electromagnetic dipole operator, 
O7, resummed to next-to- next-to-leading-logarithmic order 0, @|, and including the full 
NNLO corrections 0, S]. At lowest order, Wrilk) = 6{k). The Fm{k) are l/rrib suppressed 
subleading shape functions. In a fit to S — )■ Xg'y data only, they can be absorbed into F{k) 
at lowest order in as- The terms proportional to Cj^y are due to IR- finite bremsstrahlung 
corrections from operators other than O7. They are included at next-to-leading order (NLO) 
for 2 = 1,2,8 using the SM values for Ci^2,8- They have almost no effect on the fit, because 
they are very small in the experimentally accessible region of the photon energy spectrum. 
The coefficient C™'^' multiplying the dominant 77 contribution in Eq. (jl]) is defined as 

= Cf{^^,P^ + X:^'(^o) Q(/io) + r8(/io) Cf{^^,P^ + ■■■. (5) 
irib ^ iTib 
1=1 

Here, Cf^^iio) are the standard scheme-independent effective Wilson coefficients and rnb{fJ.o) 
is the MS 6-quark mass. The coefficients ri_6 8(/^o) contain all virtual contributions from the 
operators Oi_6,8 that generate the same effective 6 — )■ S7 vertex as O7. The ellipses denote 
included terms proportional to ln(/io/mfe) that are required to cancel the fio dependence on 
the right-hand side and vanish at fio = "rrib, such that C^"'^' is by definition /iQ-independent 
to the order one is working at. 

Since the terms in the last line in Eq. (jl]) are small, we can consider |C™^^ Vtby^*! as the 
parameter that parametrizes the normalization of the B — )■ Xg'y rate. It is extracted simul- 
taneously with F{k) from our fit to the measured E^ spectra. The important contributions 
from Oi_6.8 are the virtual corrections contained in C^"'^', which have a sizable effect on the 
normalization of the B — )■ X^'y rate. By including them in C™'^', they explicitly do not affect 
the shape of the spectrum, and so do not enter in our fit. They instead enter in the SM 
prediction for C™'^', which can be computed independently. Below, we compare to the NLO 
SM value, C™*^^ = 0.354^q;q}2 01 • For a more stringent test for new physics, evaluating C™*^^ 
in the SM at NNLO along the lines of Refs. [2I, [3| would be very valuable. 



3. FIT TO B ^ X^j 
A. Fit Setup 

To fit to the experimentally measured photon energy spectra, we insert the expansion 
for F{k) in Eq. into Eq. (jl]) and integrate over the appropriate range of E^ for each 
experimental bin and each combination of basis functions fm{x)fn{x)- The theory prediction 
for the zth bin, i?*, is then given by, 

N 

= ml \Cr%bV:,? Yl ^rnCnBi^ + ■■■, (6) 
m,n=0 

where the ellipses denote the additional included terms arising from the last line in Eq. @. 
The overall ml is expressed in terms of the moments of F[k), so it is effectively a function 
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FIG. 1: Default fit to tfie Belle and BaBar B — > Xs"f photon energy spectra from Refs. [s-ll] 
The points with error bars show the data, and the histograms show the fit result. 



of the Cn- We then perform a fit to all available bins with Co,i,...,Ar and jC'^'^'VJfeV'j*! as the 
fit parameters. We enforce the constraint Cq + ■ ■ ■ + = 1 to ensure that F{k) is properly 
normalized to unity. (An estimate of the truncation uncertainty can then be obtained from 
the size of the last coefficient.) 

As experimental inputs we use the Belle measurement from Ref. [9], and the two BaBar 
measurements from Refs. lH . The experimental statistical and systematic uncertainties 
and correlations are fully included in our fit. The BaBar spectra are measured in the B rest 
frame and are corrected for efficiencies. The experimental resolution in for each spectrum 
is smaller than its respective bin size, so we can directly use both spectra in the fit. The 
Belle spectrum from Ref. j9| is measured in the T(4S') frame and affected by both efficiency 
and resolution. Correcting the spectrum to the B rest frame depends on the shape function. 
We therefore boost our theory predictions to the T(4S') frame. Since the unfolded spectrum 
has very large bin-by-bin correlations, we apply the experimental detector response matrix 
to our theory predictions and fit to the measured spectrum. We have extensively tested our 
fitting procedure using pseudo-experiments. 

The matrices in Eq. ([H]) also have theoretical uncertainties, e.g. from higher-order per- 
turbative corrections. The corresponding theory uncertainties in the fit results are roughly 
of the same size as the experimental ones. They are not yet included in the results below. 



B. Results 



The results shown here are equivalent to those in Ref. [12|. For our default fit we use 
A = 0.5 GeV as basis parameter and four basis coefficients Cq, 1,2.3- The fit has a x^/dof = 
27.67/38 and describes the measured spectra very well, as seen in Fig. [H The fit results for 
the shape function using each of + 1 = 2, 3, 4, and 5 basis coefficients are shown in the 
left panel of Fig. [2j The corresponding results for IC^'^'^Vt^VJ*! and , where the latter is 
computed from the moments of the fitted F{k), are shown in the right panel of Fig. [2J Our 
default fit yields 

\C'r'ytbV*\ = (15.00 ± 0.54[exp]) X 10-3 . (7) 

This result is still preliminary and does not yet include theory uncertainties. It agrees within 
one standard deviation with the NLO SM value, for which we use iVt^V^*! = 40.68lo;5 ^ 10-^. 

The results in Fig. [2] verify the convergence of the basis expansion as the number of 
basis functions is increased. As one expects, the uncertainties returned by the fit increase 
with more coefficients due to the larger number of degrees of freedom. However, with too 
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FIG. 2: Fit results for two (co,i), three (00,1,2), four (co,i,2,3)) and five (00,1,2,3,4) basis coefficients and 
basis parameter A = 0.5 GeV. Left: Tlie extracted F{k) (with absorbed l/rrib corrections), where 
the colored envelopes are determined by the uncertainties and correlations of the fitted coefficients 
Cn- Right: The extracted values of |C™*^^ VtbV^*| and ml^ , where the colored ellipses show the 
respective Ax^ = 1 contours, and the gray band shows the NLO SM value for |C™'^' ^fe^tsl- 
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FIG. 3: Comparison of the fit results using different basis parameters using only two basis coeffi- 
cients (left) and five basis coefficients (right). 



few coefficients we have to add the truncation uncertainty. A reliable value for the final 
uncertainty is provided by the fitted uncertainty when the central values have converged 
and the respective last coefficients, here C3 or C4, are compatible with zero. At this point, 
the truncation uncertainty can be neglected compared to the fit uncertainties. Equivalently, 
the increase in the fit uncertainties from including the last coefficient that is compatible with 
zero effectively takes into account the truncation uncertainty. Using a fixed model function 
and ffiting one or two model parameters would thus underestimate the true uncertainties in 
the shape function. This is also seen in Fig. [31 which shows the results for different basis 
parameters A. The left plot shows the results using only two basis coefficients in the fits. 
The three fits all have a good x^/dof < 1, but disagree with each other. This shows that 
there is an underestimated uncertainty due to the basis (i.e. shape) dependence when ffiting 
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too few basis coefficients. The right plot shows the corresponding results using five basis 
coefficients in the fit. In this case, the results agree very well within the fit uncertainties. 



4. CONCLUSIONS AND OUTLOOK 

We presented preliminary results from a global fit to 5 — Xg'j data, which determines the 
total B — )• Xs7 rate, parametrized by ICy^^^VthV^*!, and the S-meson shape function within 
a mo del- independent framework. The value of jC^'^^Viby^*! extracted from data agrees with 
the SM prediction within uncertainties. From the moments of the extracted shape function 
we determine ml^ . In the future, information on rrih from other independent determinations 
can be included by a constraint on the shape function. The shape function extracted from 
B — )■ is an essential input to the determination of \Vub\ from inclusive B — )■ X^iu decays. 

A combined fit to 5 —t- Xg'j and B — )■ Xuiiy data within our framework is in progress. 
It will allow for a simultaneous determination of IC^'^^Vjby^*! and \Vub\ along with the shape 
function with reliable uncertainties. In addition to a few branching fractions with fixed 
cuts, it is important to have measurements of the B — )■ X„^z/ differential spectra (including 
correlations), e.g. the lepton energy or hadronic invariant mass spectra. As for B — t- Xg'y, 
fitting the differential spectra allows making maximal use of the B — )■ Xuiiy data, by letting 
them constrain the nonperturbative inputs and further reduce the associated uncertainties. 
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